Non-Poisson dichotomous noise: higher-order correlation functions and aging 
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We study a two-state symmetric noise, with a given waiting time distribution iP(t), and focus 
our attention on the connection between the four-time and the two-time correlation functions. The 
transition of i/'(r) from the exponential to the non-exponential condition yields the breakdown of the 
usual factorization condition of high-order correlation functions, as well as the birth of aging effects. 
We discuss the subtle connections between these two properties, and establish the condition that 
the Liouville-like approach has to satisfy in order to produce a correct description of the resulting 
diffusion process. 



PACS numbers: 02.50.Ey, 05.20.-y, 05.40.Fb, 05.70.-a 

I. INTRODUCTION 

Dichotomous noise is one of the fundamental repre- 
sentations of stochastic processes. It is used in random 
walks, quantum two-state systems, as well as other math- 
ematical models of physical and biological processes. 
This representation is used because it is simple enough to 
obtain analytic solutions to dynamical equations, yet rich 
enough to model a variety of complex physical and biolog- 
ical phenomena. The history of such two-state stochastic 
processes dates back more than a century to Markov rep- 
resentations of random telegraphic signals and yet such 
noise still finds application in models of contemporary 
complex phenomena. A few recent examples of complex 
phenomena modeled by dichotomous stochastic processes 
are disorder-induced spatial patterns ; first-passage 
and thermally activated escape Q processes; hypersensi- 
tive transport jjj; rocking rachets jah intermittent flu- 
orescence 0; stochastic resonance 0, H, Q; quan tum 
multifractality ; and blinking quantum dots l lll Il2| | . 
These and many other applications study the physical ef- 
fects of dichotomous fluctuations, either Poisson or non- 
Poisson, without addressing, however, the consequences 
that relaxing the Poisson assumption might have on the 
high-order correlation functions. 

In this paper we are interested in the high-order corre- 
lation properties of the dichotomous noise S,{t), that is, a 
symmetrical two-state statistical process with the values 
-f W and -W. Usually, for the purpose of making statisti- 
cal calculations we focus on stationary noise and use the 
stationary correlation function. 
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of realizations of the dichotomous noise. It is worth il- 
lustrating the difference between this dichotomous noise 
and a Gaussian noise with the same two-point correla- 
tion function. The difference between the two processes 
resides in the high-order correlation functions. Futher- 
more, because the noise is symmetric we only need to 
focus on even-time correlation functions. According to 
Ref. [fH , for Gaussian noise the fourth-order correlation 
function is related to the second-order correlation func- 
tion via the following expression: 

The higher-order correlation functions are analogously 
defined. In the case where all times are identical, the 
definition yields 



(e^") = (2n-l)!!(e2)". 



(3) 



a property ensuring that the distribution of ^ is a Gaus- 
sian function. By the same token, it seems natural to 
factor the fourth-order correlation function for the di- 
chotomous symmetric noise as 

{ati)at2nt3)ah)) = m)at2))mnh)), (4) 

with analogous prescriptions for the higher-order correla- 
tion functions. In the case of equal times, the definition 
Q reduces to 



(5) 



where the brackets denote an average over an ensemble 



which is similar to, but not identical to © • Equation © 
is implied for the moments of a stochastic process with 
the equilibrium distribution function 

Pi0^lm-W) + 6iC + W)]. (6) 
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Hereafter, we refer to property and the factorization 
of the corresponding higher-order correlation equations, 
as Dichotomic Factorization (DF). 

The vast majority of papers deahng with dichotomous 
noise assume the statistics of the two-states to be Pois- 
son, that is, the length of time the system remains in a 
given state has a exponential distribution. It is impor- 
tant to remark that the simplest physical phenomenon 
modeled by the stochastic variable ^ {t) is diffusion. This 
means that all the properties of the phenomenon can be 
determined by the solution to the stochastic equation 



(7) 



Allegrini et al. [Tj| found that the evolution of the 
probability density, corresponding to the dichotomous 
Langevin equation |(7|), is given by the Generalized Dif- 
fusion Equation (GDE) 



dpjx, t) 
dt 
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where the two-point correlation function under the inte- 
gral is arbitrary. 

It is interesting to note that the same GDE emerges 
from the analysis of Caceres '^^'^ studied the 

Langevin equation 



dx 



--fx{t)+at), 
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with ^(t) being a dichotomous noise and 7 a friction pa- 
rameter of arbitrary intensity. This same equation was 
studied in an earlier paper by Annunziato et al. 
It is evident that with 7 = Eq. © becomes equiv- 
alent to Eg. 0. The equation for densities found by 
Caceres [l5j is identical to that found by Annunziato et 
al. and both results for 7 ^ reduce to Eq. ((HJ. These 
results are valid independently of the form of the cor- 
relation function ^^{t). The fact that GDE is obtained 
using these different approaches is significant, since the 
work by Caceres rests on van Kampen's lemma [T^ and 
the Bourret-Frisch-Pouquet theorem J^, while the the- 
ory adopted by Annunziato et al. is the same as that 
used by Allegrini et al. Q, the Zwanzig's projection 
method jl^- In any event, both approaches adopt of a 
Liouville-like perspective. 

Bologna et al. |23| established that the GDE produces 
the same higher-order a;-moments as those derived from 
the integration of the diffusion equation, supplemented 
with the assumption that the correlation functions of 
the dichotomous variable £,{t) fit the prescription of DF. 
Bologna et al. also established that the exact solution 
of the GDE does not lead to the process of Levy diffu- 
sion, a result previously obtained using stochastic tra- 
jectories, thereby suggesting a possible conflict between 
the adoption of stochastic trajectories obeying renewal 
theory in the continuous time random walk (CTRW) for- 
malism and the adoption of a Liouville-like approach to 



the dynamics |20|. The DF assumption is not explictly 
made by Caceres However, the analysis of Bologna 
et al. indicate that the theory of Caceres implies 
the DF property. Others have also assumed non-Poisson 
statistics, while still retaining the DF property j2H. 

We establish herein that the DF condition breaks down 
as a consequence of the non-Poisson condition. Further- 
more, we show that the violation of the DF condition 
emerges from non-Poisson statistics in the same way as 
do aging properties. These results have the desirable ef- 
fect of establishing the limits of validity of the elegant 
GDE, leaving aside for the present the analysis of the is- 
sue as to whether the density and Liouville-like formalism 
are compatible with the emergence of these properties. 



II. FOUR-TIME CORRELATION FUNCTION 

In this section we show that in the non-Poisson case, 
the four-time correlation function of the dichotomous 
noise departs from the DF prescription. It has to be 
pointed out that our arguments are based on examining 
a single sequence ^, and thus on time averages, rather 
than on ensemble averages. We assume that the theoret- 
ical sequence is built up by creating a sequence {r^} of 
real positive numbers using the probability density 



(10) 



The choice of this analytical form is determined by sim- 
plicity, in which we obtain in the time asymptotic limit 
an inverse power law with index /it, while satisfying the 
normalization condition 



ll}{T)dT — 1. 



(11) 



The parameter T > insures the normalization condi- 
tion, required by the fact that ^^{t) is a probability den- 
sity and is related to the average time interval generated 
by the density. To generate a realization of the time se- 
ries we split the time axis into many time intervals of 
lengths detemined by the set of numbers {r^}. The first 
interval begins at time t — Q and ends at t — ti, the 
second begins at t = ti and ends ai t — ti + T2, the 
third begins at t = ti + T2 and ends at t — ti + T2 + T3, 
and so on. We refer to this sequence of time intervals, 
which is not observable, as the theoretical sequence. The 
dichotomic sequence under study in this paper, which 
can be observed, is created as follows. At the beginning 
of any time interval we toss a coin, and fill the interval 
with either the value W or the value —W, according to 
whether we get a head or a tail. Thus, if we move along 
the observable sequence, we meet large time portions of 
the sequence, within which the sequence retains the same 
value, either W or —W. We refer to these time intervals 
with the same value of ^, as experimental laminar regions 
and to the corresponding distributions of time lengths as 
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ipexpiT). The adoption of the suggestive term experimen- 
tal reflects the fact that this procedure is the same as the 
one we would adopt when making a real experimental 
observation. A relevant example is the phenomenon of 
blinking quantum dots |12j ■ which has been the object 
of some very interesting theoretical papers [2^ using 
dichotomous stochastic processes. A single quantum dot 
undergoing a process of resonant fluorescence produces 
an intermittent light signal, which can be identified with 
the sequence ^(t) here under study, with W and —W, 
meaning light-on and light-off, respectively. 

We point out that 'ipexpi'r) does not necessarily coincide 
with'(/;(r). According to Ref. [2^ the theoretical waiting 
time distribution ip{t) is connected to the experimental 
waiting time distribution by the Laplace transform rela- 
tion 

^(.) = ^M^, (12) 

where the Laplace transform of a function / {t) is denoted 
by / (u). However, in the time asymptotic limit ipexpi'T') 
has the same inverse power law form as does ipir), that 
being Eq. (|1(J|I . with the same power-law index /i. In the 
special case of blinking quantum dots the experimental 
waiting time distribution is found to be an inverse power 
law with index fi < 2. Here we consider the complemen- 
tary case fj, > 2, so as to realize a condition compatible 
with the existence of a stationary correlation function for 

Due to the theoretical prescription that we adopt to 
realize the dichotomic sequence under study, a given ex- 
perimental laminar region, namely, a time interval where, 
as earlier pointed out, ^(t) keeps the same sign, might 
correspond to an arbitrarily large number of theoretical 
time intervals, to which the coin tossing procedure as- 
signs the same sign. We shall refer to these theoretical 
time intervals as theoretical laminar regions, or, more 
simply, as laminar regions. It is evident that the begin- 
ning of a laminar region corresponds to the occurrence 
of a random event, namely the coin tossing that deter- 
mines its sign. The laminar regions are not observable, 
while the experimental laminar regions are observable, 
by definition, and begin and end with a random event. 
We cannot establish if other random events occur or not, 
and how many, between the beginning and the end of an 
experimental laminar region. 

The theoretical approach that we adopt in this section 
rests on the same time average procedure as that adopted 
by Geisel et al. . Let us devote some attention to the 
prescription given by these authors to evaluate the two- 
point correlation function $^(1^2 — ^il) psj : 

oc 

/ [T-{h-t,mT)dT 

'^dh-h)^ ''-'' , (13) 

Jo Ti;(T)dT 

where we assume ^2 > ^i- This equation for the cor- 
relation function implies that, with a window of size 



A = <2 — ii we move along the entire (infinite) theo- 
retical sequence of laminar regions and count how many 
window positions are compatible with the window being 
located within a theoretical laminar region, which must 
have a length larger than the window size. In addition we 
have to count the total number of window positions. In 
other words, the stationary correlation function of ^ (t) is 
nothing but the probability that the two times ti and t2 
are located within the same laminar region. If these two 
times are located in different laminar regions, the adop- 
tion of the coin tossing procedure for any contribution of 
a given sign to the correlation function would produce, 
with equal probability, a contribution with opposite sign, 
thereby providing a vanishing contribution. An attrac- 
tive way to explain this procedure is through the concept 
of random events. First of all, the lengths of the lami- 
nar regions are determined by the random drawing of 
the numbers r, with distribution 'iP{t). At the border 
between one laminar region and the next we toss a coin 
to decide the sign of the next laminar region. This coin 
tossing is a random event and no random event can occur 
between two times located in the same laminar region. If 
the two times are located in different laminar regions, 
one or more random events must have occurred between 
them. Thus the correlation function 'I'^dti — ^2!) can also 
be interpreted as the probability that no random event 
occurs between times ti and ^2- 

We evaluate the four-time correlation function, using 
the same arguments. Consider four times, ordered as 
< t2 < is < t^. The corresponding correlation func- 
tion exists, under the following conditions. The first con- 
dition is that all four times are located in the same lami- 
nar region. The second condition is compatible with the 
pairs (ti, and (ta, t^) being located in distinct laminar 
regions. This means that the times ti and ^2 belong to a 
laminar region, denoted by Ti^2, the times ^3 and be- 
long to a laminar region denoted by 23,4, and Ti_2 7^ 73,4. 
Using the random event concept, the second condition 
implies that no random event occurs between ti and i2, 
or between ^3 and t4, while at least one random event 
occurs between t2 and ^3. 

We use the notation p{ij) to denote the probability 
that ti and tj belong to the same laminar region. Thus 
the prescription for the correlation function given by Eq. 
(|13|l can be expressed as the probability function 

$?(^2-ti)=p(12). (14) 

We also use the notation 

p(i]) = l-p {ij) (15) 

to denote the probability that at least one transition oc- 
curs between times ti and tj. It is convenient to use the 
conditional probability concept, and the Bayesian nota- 
tion (see, for instance, H^). We denote the joint proba- 
bility of events A and B by p{A, B) and the conditional 
probability of occurrence of event A given event B with 
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p{A\B). Thus, we have 

p{A\B) 



p{A,B) 
p{B) 



(16) 



We denote the conditional probabihty that event A oc- 
curs, given that event B does not, by p{A\B). Using the 
prescription of Eq. , the latter conditional probabil- 
ity, p{A\B), is expressed as follows 



p{A\B) 



p{A)~p{A,B) 
l^p{B) ■' 



(17) 



where we have used the relation p {A) = p(A, iJ)+p(A, B) 
for the numerator. 

The probability that times ti and tj belong to the same 
laminar region j and that, simultaneously, times tr and 
is belong to the same laminar region T^ g^ regardless of 
whether T^j coincides with T^^s, or not, is a joint prob- 
ability expressed by the symbol p(ij, rs). Thus the four- 
time correlation function can be formally expressed as 
follows: 



(e(ii)e(i2)e(i3)e(i4)) 



2\2 



:p(12,34). 



(18) 



On the other hand, using the notation introduced ear- 
lier, we have two contributions to the four-time correla- 
tion function. The first contributions is determined by 
all four times being in the same laminar region with no 
random event occuring between t\ and (condition 1), 
whereas the second contribution corresponds to the prob- 
ability that at least one random event occurs between t2 
and ^3, given the condition that no random event occurs 
between t\ and ti and none between t-^ and t/^ (condition 
2): 



p(12,34) =p(14) -|-p(23)p(12|23)p(34|23). 



(19) 



Eq. (|19|l corresponds to the superposition of independent 
contributions from condition 1 and condition 2. The con- 
tribution due to condition 1, p(14), according to the ear- 
lier definitions, is the probability that t\ and ^4 belong to 
the same laminar region. The contribution due to con- 
dition 2 is given by the second term on the right hand 
side of Eq. I|19|l . Again, according to the notation that 
we are using, see Ea. (ll5f) . p(23) is the probability that a 
random event occurs between and ^3, thereby discon- 
necting the two laminar regions. Consequently p(12|23) 
is the probability that t\ and belong to the same lam- 
inar region given that at least one random event occurs 
between and t3. Finally, p(34|23) is the probability 
that is and t4 belong to the same laminar region, given 
that at least one random event occurs between ti and 
t-i- Thus, the product of these three probabilities is the 
appropriate quantity corresponding to condition 2. 

To transform the equality Eq. 1)191) into a relation 
involving correlation functions, we use Ea. (|18|l . for the 
four-time correlation function. The two-time correlation 
functions emerge from the second term on the right hand 



side of Eq. ((T^ via the proper use of Eq. ((T^ , Eq. lfT7|) 
and Eq. H15|l . Thus, we obtain 



$e(i4-ii) 



(20) 



($4fe - ii) - ^cfe - <i))(^cft4 - H) - - ti)) 

l-$^(i3-t2) 

Eq. (|20|l is a major result, being an exact expression for 
the four-time correlation function independently of the 
statistics of the dichotomous process. We stress that the 
general form of Eq. (|2()(l is not factorable and is therefore 
distinct from DF. 

Note that in the Poisson case, the waiting time distri- 
bution tpit) is exponential. Using the prescription given 
by Eq. 1)13(1 it is not difhcult to show that the correla- 
tion function of ^ is also exponential. Then, after tedious 
but straightforward algebra, we establish that Eq. H2U|) 
reduces to 

{ati)at2)at3)au)) = {ati)at2)){ah)au)), (21) 

which coincides with Eq.Q), that is, the process becomes 
compatible with the DF. Given that the DF holds true 
for the four-time correlation function, it is possible to ex- 
tend the DF property to the 2A^-time correlation function 
using induction. 

Thus, we conclude that the four-time correlation condi- 
tion 1)21) [I . for waiting times that have non-Poisson statis- 
tics, violates the DF underlying Eq. ©. This violation of 
the factorization property seems to be a satisfactory ex- 
planation of why the GDE pO] does not yield the proper 
Levy diffusion in the asympotic limit. On the other hand, 
using the results of this section we recover the results of 
the numerical calculations and theoretical prediction of 
the fourth moments obtained by Allegrini et al. |23| . To 
establish this latter point we integrate Eq. Q with the 
initial condition x{Q) for all the trajectories. Further- 
more, we evaluate the fourth power of x{t), and average 
over all the trajectories of the Gibbs ensemble. By us- 
ing the stationary condition, which makes this correlation 
function depend only on the time differences, rather than 
on the absolute time, we obtain 

{x\t)) (22) 

t t4 t-i ±2 

dU J dh J dti j dti{^{ti)^{t2)at3)^{t4)). 


Introducing the newly obtained expression for the fourth 
order correlation Eq. I)2UI) into 122)l . in the time asymp- 
totic limit the leading contribution to the fourth moment 
is given by the first term on the right hand side of Eq. 
()20)l . Therefore we replace the integrand in Eq. 122)) 
with <i>^(i4 — ii), and using the inverse power-law form 
of the correlation function, we carry out the four time 
integrations and obtain {x'^{t)) cx i^"''. By extending 
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this way of proceeding to the calculation of the 2n-times 
correlation function, we derive the general result 

cxi2"-^+2, (23) 

for 2 < ^ < 3, and {x'^'^{t)) oc t^"^^ for ^ > 3, in agree- 
ment with the numerical results of Ref . 27] . 

The asymptotic result H23|l establishes that the 2n- 
moments do not have the scaling corresponding to the 
DF condition. If we assume that the condition of Eq. 
Q applies, in keeping with the nature of the GDE, in- 
stead of (1231 we would obtain (a;2"(t)) oc t2n(4-M)/2^ ^j^j^ 
one factor of /i occuring for each order of the moment. 
Consequently, the DF implies the existence of the scaling 
X (xt^ , with the scaling index given by 

6 = for 2 < /.t < 3, J = i for ^ > 3, (24) 

where /i — 1 is the Levy index. This later result agrees 
with the scaling predicted by the GDE, as established in 
Ref. 20]. Here the central fact to keep in mind is that 
Eq. (fTI) generates Levy walks, rather than Levy flights. 
A Levy flight is a kind of random walk in which the step 
lengths have an inverse power-law distribution, so the 
second moment of the dynamical variable diverges. The 
Levy walk, on the other hand, ties the length of a step 
to the time required to take the step, resulting in a finite 
second moment for the dynamical variable. Furthermore, 
it takes an infinite time for a Levy walk to yield the same 
scaling as a corresponding Levy flight, the latter scaling 
index being given by 

5 = for 2 < /i < 3, (5 = i for ^ > 3. (25) 

/i — 1 2 

For this reason, the Levy walk, introduced by Shlesinger 
et al. [25 , can be considered to be a manifestation of the 
Living State of MatteWLSM) |2^ , in the sense described 
in some recent work j3fi l3lj . The LSM is interpreted as 
the existence of a scaling condition intermediate between 
that of dynamics and thermodynamics and which can 
last forever. 



III. AGING 

In this section we adopt the Bayesan formalism to eval- 
uate the correlation functions in a non-stationary condi- 
tion. This enables us to establish that the breakdown of 
the DF condition is closely related to aging. 

Before proceeding with the formalism, we briefly re- 
view why non-Poisson statistics produces aging, as dis- 
cussed in detail in Refs. ji2j. Suppose that we create 
an infinite sequence of time intervals of length , namely, 
the theoretical sequence discussed earlier. As mentioned, 
we create the observable sequence by filling the time in- 
tervals, called laminar regions, with either W or —W, 
according to the coin tossing prescription, with the first 



laminar region beginning at time t = t^. Let us imag- 
ine, to facilitate the discussion of this section, that the 
theoretical sequence is observable, even if in practice it 
is not. If we begin the observation process at the same 
time when the theoretical sequence is generated, the re- 
sult of our observation yields the waiting time distribu- 
tion of Eq. (|10|l . If the observation of the theoretical 
sequence begins at a given time ti > tg, the distribution 
of the waiting times before the first exit from the laminar 
region, denoted by ipti,toit), will not coincide with ipit)- 
This is a consequence of the first laminar region observed 
having begun at any time between ti and Iq. Thus, the 
resulting waiting time will be, in general, shorter than 
the real sojourn time generated by iP{t). In the Pois- 
son case this shortening of the time does not have any 
effect on the shape of ^Jti,to{t), which remains identical 
to iP{t). In the non-Poisson case, on the contrary, delay- 
ing the process of observation does influence the shape 
of .tn (t) causing it to depart from the form of i^ir) 

Let us now address the problem of building up the ag- 
ing correlation function of ^{t). We study the correlation 
between ^(t2) and £,{ti), with the condition that ^2 > > 
to] to being the time at which the laminar region begins. 
We solve this problem in two steps. In the first step we 
define the correlation function A^*°\t2 ~ ti), without re- 
quiring that the laminar region begins a,t t = to, but that 
it in fact begins at a time intermediate between ti and 
to- This corresponds to stating that A'^*°^(t2 — ii) is a 
correlation function of undefined age, younger, though, 
than the (ti — to)-o\d correlation function. In the second 
step we set the additional condition that the laminar re- 
gions begin a,t t = to, and we give the prescription to 
determine the correlation function a notation de- 

noting in fact the (^i — to)-old correlation function. The 
latter aging correlation function fits the earlier definition 
of "0*1,40 (^)- The corresponding analytical expression will 
make it possible to establish the effect of aging on the 
phenomenon, namely the effect of moving both t2 and 
ti away from to as well as the more traditional effect of 
increasing the distance between t2 and ti. 

Note that the first correlation function is given by 

A(*«'(t2-ii)=p(A|B). (26) 

This identification of A'-'^°^{t2 — ti) is consistent because 
we define A by the condition that both ti and t2 belong 
to the same laminar region, while B is defined by the con- 
dition that ^0 does not belong to the same laminar region 
as ti. Of course, with this interpretation B is defined by 
the condition that to belongs to the same laminar region 
as ti. 

The conventional correlation function is the probability 
that ti and t2 belong to the same laminar region, and 
thus is the probabilty that property A occurs, so we can 
write the second equality 

$?(t2~ii)=p(A). (27) 
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The probability that A and B take place in the same lam- 
inar region allows us to write, in terms of the correlation 
function, 



^dt2 - to) ^ piA, B) 



(28) 



In fact, this is the probability that t2 and Iq belong to 
the same laminar region, and, thanks to the time order- 
ing ^2 > ii > ^Oi this is equal to the probability that 
both A and B occur. Finally, the probability that ii and 
to belong to the same laminar region, namely, the prob- 
ability that the property B applies, enables us to write 



*e(ii -<o) =p(5)- 



(29) 



At this stage, to express A^*°)(<2 — ti) in terms of more 
familiar correlation functions we insert Eq. H17|l into 
Eq.ljSnj, which yields 

i - 9^(ti - to) 

It is easy to show that in the Poisson case Eq. H3U|) reduces 
to 



A(*°)(i2-ii) = $e(t2-ti), 



(31) 



independently of to- 

Now let us take the second step, and explicitly evaluate 
^f''\t2-h). This aging correlation function is the sum 
of two probabilities. The first contribution is the prob- 
abiliy that no event occurs between to and t2, thereby 
ensuring that ti and t2 belong to the same laminar re- 
gion. The second contribution is the probability that an 
arbitrary number of events occurred between to and t2. 
Note that the laminar region beginning at t = to implies 
that at this time a random event occurs, which is in fact, 
the beginning of the laminar region. As stated a number 
of time earlier, at the beginning of any laminar region, 
we toss a coin, to decide the sign of the laminar region. 
This is the random event that makes it possible for us to 
express *f*°^(t2 -ti) as follows 



<i>f''^(t2-ti) = *(t2-to) 



^^^£(^2 - tl) - $f(t2 - to) , , 

+ l-^,(tr-to) 

In Eq. (|32() we have used the conventional notation of 
the CTRW formahsm 



dtV(t'), 



(33) 



where ipit) is the waiting time distribution of Ea. (|10|l . 
MontroU and Weiss [s^ make the implicit assumption 



that the laminar region begins at t = 0. Thus, \E'(t) is 
the probability that no event occurs up to time t, after 
the random event occurs at time t = 0. Here we replace 
the initiation time t — with t = tg. Thus, ^'(t2 — to) 
is the probability that no random event occurs between 
to and t2, as required. The second term in Ea. (|32(l is 
the product of the probability that one or more events 
occurred between ti and to, given the fact that t2 and ti 
are in the same laminar region and to is not. 

We note that Eq. 1)32(1 interrelates factorability and 
aging and consequently is the most relevant expression 
for our discussion. The importance of this result can be 
made transparent by going back to the discussion in Sec- 
tion 2. Eq. (|20|l . the expression for the fourth-order cor- 
relation function, can be reexpressed as, using Ea. (|30|l . 



= $^(t4-ti) 



+ ($4(t2 - ti) - $^.(t3 - ti))A(*^)(t4 - h). (34) 

As pointed earlier, in the Poisson case, see Eq. (|31|l . 

A^''\U~t3) = <P^{U~t3), (35) 

independently of t2. By inserting Ea. H35|l into Eq. H34|l. 
and noting that ^^{t4 - ti) <i>^(t4 - t3)<i>^(t3 - ti), we 
see immediately that the DF condition is recovered: 



= 4'«(t2 -ti)$5(t4 -ta). (36) 



Thus, we have established that the breakdown of the DF 
condition and aging are interrelated. In fact, annihilating 
the aging property has the effect of reestablishing the DF 
property. 



IV. CONCLUDING REMARK 

The equivalence between the trajectory and density 
pictures of physical phenomena is one of the major ten- 
ants of modern physics. It therefore came as quite a 
suprise when Bologna et al. discovered an incon- 

sistency between these two pictures in the case of non- 
ordinary statistical mechanics. The form of the inconsis- 
tency had to do with the derivation of anomalous diffu- 
sion of the Levy kind, using dichotomous noise and either 
CTRW or the generalized central limit theorem. Both of 
these approaches use trajectories and not the Liouville- 
like approach for densities, such as does GDE. It is a 
simple matter, using Eq.©, to show that GDE yields a 
hierarchy of moments (a;^" (t)) with n — 1,2, which 
coincides exactly with the hierarchy generated by fluctua- 
tions ^ (t) satisfying DF. This factorization, obtained us- 
ing the density, contradicts the hierarchy generated using 
the trajectories in Section 3. We have limited the anal- 
ysis to the fourth-order correlation functions, however. 
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this order is sufficient to identify the source of the incon- 
sistency between the trajectory and density pictures as 
being due to the non-Poisson character of the statistics. 

We have also shown that a departure from Poisson 
statistics has the effect of introducing a memory into 
the correlation functions that can last for an infinitely 
long time. For dichotomous noise the two-time corre- 
lation function, using either trajectories or densities is 
the same, however, higher-order correlations are not the 
same for non-Poisson statistics. The deviation from Pois- 
son statistics is manifest in a dependence of correlations 
on the difference between the initiation time and the ob- 
servation time, that is, on the age of the system. Age 
destroys the DF property and may represent a state of 
matter intermediate between the dynamic and thermody- 
namic condition, mentioned earlier, the Living State of 
Matter. This eternal state of nonequilibrium, in which a 
perturbed phenomena relaxes to, but never attains, equi- 
Ibrium, should be contrasted with the Onsager Principle 
in which physical systems are assumed to be aged. An 
aged physical system is one that has reached equilibrium 
with a heat bath long before measurements are taken. 

It is evident that to establish a density picture equiv- 
alent to the trajectory picture, in which the time aver- 
ages and ensemble averages are the same, in the non- 
Poisson as well as in the Poisson case, we have to over- 
come the limitations of the Liouville-like approaches of 
Refs. [13 El El- This difficult issue calls for further re- 
search. Nevertheless, the merit of the present paper lies 
in the fact that it has revealed the violation of the DF 
property when the statistics of the underlying process are 
non-Poisson. DF is a factorization property assumed for 
dichotomous noise by researchers in multiple fields, often 
without the realization that such factorization is tied to 
the statistics of the process. 

It is worth remarking that Eq. the general diffu- 
sion equation, can also be derived by assuming that the 
random sequence ^ (t) is built up by time-modulating a 
generating Poisson distribution ^{X{t),t) — exp{—X{t)t) 



as shown in detail by Bologna et al. [s^. The result- 
ing sequence, however, is not a renewal sequence, such as 
found in CTRW. We need to understand why abandoning 
the Poisson assumption and adopting a Liouville-like ap- 
proach leads to physical effects that are inconsistent with 
renewal processes. This is a difficult problem whose so- 
lution also requires additional research. It is important 
to point out, to avoid any possible confusion, that the 
GDE is widely used to describe transport processes (see 
Refs. H^H^, for some recent papers). However, these 
papers refer to subdiffusion, a physical condition where 
the correlation function of the fluctuation f cannot be 
defined; not even in the non-stationary sense of Section 
3. The discussion herein focuses on superdiffision and 
addresses the problem of computing high-order correla- 
tions for renewal process with non-exponential waiting 
time distributions. The solution to this problem is given 
by Eq. (|20|l . however this crucial property has not yet 
been obtained using Liouville-like methods 17, 18, 19]. 

In conclusion, by means of the conditional probabil- 
ity formalism, we have found the exact expression for 
the fourth-order correlation function, and we have shown 
that in the non-Poisson case, this expression violates the 
DF condition. We have also established a close connec- 
tion between the DF breakdown and aging. In the case 
where ^ > 2 the aged condition is possible. However, if 
/i < 3, the aging condition lasts forever. We see, in fact, 
from Eq. H13() that in this case the correlation function 
^^(t) is an inverse power law with index /i — 2. Thus, it 
takes an infinitely long time for the age-dependent cor- 
relation function Eq. H32|l to become stationary. This 
is a remarkable result, which challenges the traditional 
treatments of such stochastic dynamical processes based 
on the generalized master equation (GME). The analy- 
sis of the GME based on these results will be taken up 
elsewhere. 
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